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ABSTRACT. We show that the evolution of the universe is sin- 
gularity free in a class of graviton-dilaton models. 
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We present a general analysis of the evolution of a homogeneous isotropic 
universe in graviton-dilaton theories and show that the evolution of the ob- 
served universe is singularity free in a class of models. Q 

1. Action and Equations of Motion 

In the absence of a dilaton potential, the most general graviton-dilaton 
action is given by 

S = -\J d A x^~g (xR + ^ (V X ) 2 ) + S M (M, 9iat ) , (1) 

where x is the dilaton with the range < x < 1 u(x) is the arbitrary 
function that characterises the theory. In ([]]), the matter fields feel only the 
gravitational force and, hence, the physical metric is g^ v and the physical 
quantities are those directly obtained from g^ u . 

In the model derived in 0, the function fl(x) = 2w(x) + 3 is required to 
satisfy 

o _ 3 

finite Vn>l, 0<x<l 
oo at x = 1 only 

^i(l-xr 2a , ^<a<!> (2) 

where Qq > and Q\ > are constants. We further assume, for the sake 
of simplicity, that is a strictly increasing function of x- The function 

Q(x) is otherwise arbitrary. 

In the following, we take "matter" to be a perfect fluid with density p and 
pressure p, related by p = "yp where —1 < 7 < 1 and take the line element 
to be given by 

rfs 2 = _ dt 2 + e 2A(t) ^2 + r 2( d6 2 + ^2 (3) 

1 The possibilties of such a model arising within string theory is discussed in [Q, [|. See 
also JH . Other proposals for singularity free evolution of the universe in string theory can 
be found in Q. 

2 \(today) — 1 corresponding to a non zero value of the Newton's constant. 
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The equations of motion from ([!]) are 
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where upper dots denote t-derivatives, e = ±1 and we have used Q = 2oj + 3 
and p = '-/p. The square roots are to be taken with a positive sign always. 
Equations @ and (JD can also be written as 



m 

dx 



-3 A 



(a(t) + c) 



-1 + e sign(x) VK , 
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where c = x\/^e 3j4 \t=ti is a constant, ti is an initial time, and 



a(t) 
K = 
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t e -^ A 
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Note that under time reversal t — > —t, e — > — e in equation whereas 
equation (RI) is unchanged. Thus, evolution for — e is same as that for e, but 
with the direction of time reversed. 

Note also that, in general, there will appear in the solutions positive 
non zero arbitrary constant factors e A °, po, and %o i n front of e A , p, and 
X respectively. However, it follows from ([|) and @ that they can all be 
set equal to 1 with no loss of generality if time t is measured in units of 
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Therefore, in the solutions below we often assume that t is 



measured in appropriate units and, hence, set these constants equal to 1. 
Note, however, that the model dependent constants associated with fl(x), 
such as Qq or Oi in (El), cannot be set equal to 1. 



2. Analysis of the Evolution 
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Our main goal is to determine whether the evolution of the universe in 
the present model is singular or not. The task is trivial if one can solve the 
equations of motion (^) - (^) for arbitrary functions (or i/j (</))). However, 
explicit solutions can be obtained only in special cases [§, |l| and the methods 
involved in obtaining such solutions are inapplicable in the general cases of 
interest here. Also, even when applicable, the details of the solution tend to 
obscure the general features of the evolution. Hence, a different approach is 
needed which is valid for any matter and for any arbitrary function £l(x) and 
which enables one to analyse the evolution for the general case and obtain 
its generic features even in the absence of explicit solutions. We will present 
such an approach below. 

Note that for the singularities to be absent, all curvature invariants must 
remain finite. A sufficient condition, proved in || , for all curvature invariants 
to be finite is that the quantities 

e- ; &,JL ; * «£ and V„>1 (11) 

X X*l X X « dx n \xj 

be all finite. To proceed, we also need initial values of A, x, and X or 
equivalently Q at initial time tj. However, their numerical values are not 
needed in our approach as long as they are generic f] which we assume to be 
the case. 

The value Q(U) £ 3 will turn out to be special. Thus, there are 16 sets of 
possible initial conditions: 2 each for the signs of Afa), x(U)) and e, and 2 
for whether fi(ij) < 3 or not. 

However, it is not necessary to analyse all of these 16 sets. The evolution 
for — e is same as that for e, but with the direction of time reversed. Hence, 
only 8 sets need to be analysed. Now, let e = 1. If % < then A can not 
be negative for any value of Q, see ([|). Similarly, if % > and Q > 3 then 
also A cannot be negative since \ff- + — y§ > m equation (^). 
Thus, there remain only 5 sets of possible initial conditions. We choose these 
sets so as to be of direct relevence to the evolution of observed universe and 
analyse each of them for t > tj. We then use these results in section 6 to 
describe the generic evolution of the universe. We also assume that p ^ 
and 7^1 identically, as is relevent to the observed universe. 

3 A non generic example: x(ti) = and xi^i) — b *- e - ®\u = °°- 
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Two remarks are now in order. First, the amount and the nature of 
the dominant "matter" in the universe are given by po an d 7 where — 1 < 
7 < 1. The later varies as the scale factor e A of the universe evolves. If e A is 
increasing, i.e. if A > 0, then the universe is expanding, eventually becoming 
dominated by "matter" with 7 < |, when such "matter" is present. For the 
observed universe, which is known to contain dust for which 7 = 0, one may 
take 7 < 0. Thus, by choosing t,i suitably we may assume, with no loss of 
generality, that 7 < | (< for observed universe) if A > 0. Similarly, if 
A < then we may assume, with no loss of generality, that 7 > | when such 
"matter" is present. This is valid for the observed universe which is known 
to contain radiation for which 7 = |. 

Second, consider a(t) in (|7]). The integrand is positive and oc -7=. (The 
dependence on p can be absorbed into the unit of time, as explained in 
section 3.) Hence, the value of the integral is controlled by Q. For example, 
in the limit x ~ * lj it is controlled by the constant Q\ in (0). We will use this 
fact in the following. Also, since t > U, the sign of a(t) is same as that of 
(1 — 37). Moreover, a(t) = if po = or 7 = | corresponding to an universe 
containing, respectively, no matter or radiation only. 

2 a. A(ti) > 0, x(U) > 0, and flfo) > 3 

These conditions imply that e = +1 in (M and the constant c > in (|7]). 
For t > ti, the scale factor e A increases and, eventually, 7 can be taken to be 
< I when such "matter" is present. In fact, 7 < for the obeserved universe. 
Then, (1 — 37) > and, hence, a(t) increases. Thus, for t > t i? \ > and x 
and Q both increase. Since Q > 3, it follows from equation (||) that A > for 
t > ti and, hence, e A increases. Eventually, x ~ > 1 an d — > — x)~ 2a - 

In the limit as x ~* 1) ( a (t) + c ) > is a constant to an excellent 
approximation. Then, equation (|8[) can be solved relating x an d e A : 

e A = (constant) (1 - x)~^^ • (12) 

Substituting this result in equation (|7]) then yields the unique solution in the 
limit x ~ > 1 ; 

e A = e A °t3(TT7) ; x = 1 _ x £~ (i+^u-oo , (13) 

where Aq and %o > are constants and t is measured in appropriate units. 
Since % — > 1, it follows from ([HJ) that t — > 00. It can now be seen that 
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(cr(£) + c) > is indeed a constant to an excellent approximation. Also, the 
quantities in (|TTD are all finite, implying that all the curvature invariants are 
finite. Thus, there is no singularity as x ~~ * 1 an d, thus, for t > t{. 

2 b. A(U) > 0, x(U) < 0, fl(tj) > 3, and e = 1 

The constant c < in (0). For t > U, x an d ^ decrease and the scale 
factor e A increases. Equation (0) now becomes 

ax 

where K(t) is given in 

Consider first the case p = or 7 = 1, discussed in section 2a and which 
led to the present phase. Then, p (l — 37) < and, therefore, a(t) < 0. 
Thus, (&(t) + c) < c < and, hence, x(£) is negative and non zero, implying 
that x(t) decreases for t > £j. Also, K(t) remains finite and it follows from 
equation ([8]) that ^ < 0. Thus, A > since x < and, hence, e A increases 
for t > U. 

Eventually, as t increases, x — > and e A — > 00. Also, 2p °*^) +e )2 <C 1 in 
(H) since p = or 7 = 1. Equation (§) can then be solved relating x an d 
e A : _ 

e = (constant) X 6 , (14) 

where x ~ y and Qo is given in (||). Substituting this result in equation (0) 
then yields the unique solution in the limit x ~* 0: 

e A = e Ao (t-t ) n , X = Xo(t-to) m , (15) 
where A , Xoj an d > U are some constants and 



3 + V3n -2 

n = — . , m = ; . (lb) 

3(1 + y/m^ 1 + 

Note that m < 0. Then, since x ~~ > 0, it follows from (|15D that i — > 00. 

Thus, as x — > 0, £ — > 00 and e A — > 00. It can also be seen that the 
quantities in (]XT|) are all finite for £; < £ < 00, implying that all the curvature 
invariants are finite. Thus, there is no singularity for ti < t < 00. 
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Consider now the case where p and 7 are non zero and have generic values. 
We have A(U) > and x(U) < 0. Hence, e A increases and x decreases. As 
e A increases, 7 can eventually be taken to be < h, when such "matter" is 
present. In fact, 7 < for the obeserved universe. In such cases where 7 < 0, 
a(t) grows faster than t as can be seen from ([|). Then, as t increases, the 
factor (a(t) + c) becomes positive. 

This implies that x, initially negative, passes through zero and becomes 
positive. We then have A > and x > 0. As t increases further, x an d, 
hence, fi(x) increase. If Q > 3 then further evolution proceeds as described in 
section 2a. If O < 3 then, upon making a time reversal, the initial conditions 
become identical to the one described in section 2e. The evolution in the 
present case then proceeds as in section 2d for the case which leads to the 
initial conditions of section 2e. 

2 c. A(U) < 0, x(U) < 0, and fi(fc) > 3 

These conditions imply that e = — 1 in (0J) and the constant c < in 
(0). For t > ti, the scale factor e A decreases and, eventually, 7 can be taken 
to be > | when such "matter" is present. In fact, this is the case for the 
observed universe. Then (1 — 37) < and a(t) decreases or remains constant. 
Therefore, the factor (cr(t) + c) < c < and, since e A decreases, we have that 
x(t) < x{U) < for t > ti. Thus x an d, hence, Q decrease for t > U. Also, 
^ > since x < and A < 0. 

In (H), K{ti) > 1 since Q > 3. From the behaviour of e A , x, fi(x), and 
(cr(t) + c) described above, it follows that K decreases monotonically. The 
lowest value of K is ^ < |, achievable when x vanishes. This, together with 
K(ti) > 1 and the monotonic behaviour of K(t), then implies that there 
exists a time, say t = t m > U where K{t m ) = 1 with x(tm) > 0. Hence, 
T^-(^m) = 0. Therefore, A{t m ) = since x(t TO ) is non zero. This is a critical 
point of e A and is a minimum. Also, equation (§) can be written as 

A{U) - A(t m ) = T (tl) p (-1 + y/K) = finite , (17) 

J X(tm) ^X 

where the last equality follows because both the integrand and the interval 
of integration are finite. Therefore, A(t m ) is finite and > —00 and, hence, 
e A{t m ) j s finite and non vanishing. 
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Thus, for t > ti, the scale factor e A continues to decrease and reaches a 
non zero minimum at t = t m . The precise values of t m , A(t m ), x{tm), an d 
Q(t m ) are model dependent. Hence, nothing further can be said about them 
except, as follows from K(t m ) = 1, that £l(t m ) ^ 3 in general and = 3 when 
p = 0. 

However, the above information suffices for our purposes. It can now be 
seen that the quantities in ([U]) are all finite, implying that all the curvature 
invariants are finite. Thus, there is no singularity for U <t < t m . 

For t > t m , one has A(t) > 0, x(t) < 0> an d ^(x(tm)) ~ 3 by continuity. 
Further evolution is analysed below. 

2 d. A(U) > 0, x(U) < 0, fl(ti) £ 3, and e = -1 

This implies that the constant c < in (0) and that ^(U) < in (|8|). 

Hence, K(U) < 1. For t > U, the scale factor e A increases and, eventually, 7 
can be taken to be < | when such "matter" is present. In fact, 7 < for the 
observed universe. Then, (1 — 37) > and, hence, a(t) increases. Thus, e A 
is increases, x decreases, and (cr(t) + c) increases. Now, depending on p , 7, 
and the details of fl{x), K(t) in (^|) may or may not remain < 1 for t > ti. 

Consider first the case where K(t) remains < 1 for t > U. It then follows 
that (cr(t) + c) cannot have a zero for t > tj and, since c < 0, must remain 
negative and non infinitesimal. Note that 7 must be > for otherwise a(t) 
grows as fast as or faster than t, making (<r(t) + c) vanish at some finite time. 
Let 

lim (a(t) + c) = c e (18) 

t — >oo 

where c e is a negative, non infinitesimal constant. We then have for <t < 
00, c < (cr(t) + c) < c e < and, hence, x{t) < 0- Therefore, eventually 
X(t) -> 0. Since K(t) < 1 and £(t) < 0, it also follows from (g) that g < 

and, hence, A(t) > 0. We will now consider the limit x ~~ > 0- 
K(t) < 1 implies that as x — ^ 0, A' ^ ^ < 1 where 

is a constant and f2 is given in @. Now, as x ~ > 0, equation (||) can be 
solved relating x an d e A : 

e A = (constant) x ~ 3_% f n " _ ( 2 0) 
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Note that (3 — \/3fie) > since Q e < 3. Hence, e A — > oo as x ~^ 0. 
Substituting (p0| ) in equation (^) then yields the unique solution in the limit 
X -> 0: For fi e ^ §, 



e = e (t - sign(m)t) , x = Xo (<o - sign(m)t) , (21) 

where t is measured in appropriate units, A , Xo, an d to > U are some 
constants, and 



n 



y/3TT e -2 



3(1 - y/3TQ ' 1 - 



m = == . (22) 



e A = e A °e- C -^, X = Xoe^K (23) 



For n e = §, 



where c e is defined in (|n_ 

Let Vt e > |. Hence, m > 0. Also, n < because fl e < 3. Since x ~ y 0, it 
follows from (|2l| ) that t — > to > U> which implies that x vanishes at a finite 
time to- In this limit, the scale factor e A — > oo since n < 0. The quantities 
in (0), for example 2£, also diverge, implying that the curvature invariants, 
including the Ricci scalar, diverge. Thus, for D, e > |, there is a singularity 
at a finite time to- 

Let f2 e < |, which is the case in our model, see (0). When f2 e < |, m < 0. 
Also, n > because f2 e < 3. Since x — ► 0, it follows from (|2l| ) that t — > oo. 
In this limit, the scale factor e^ 4 — > oo since n > now. The same result 
holds also for the solution in (|23"D with J7 e = ~. 

For -ftT(t) to be < 1, (cr(t) + c) must not be too small. Since c < 
and cr(t) > 0, it necessarily implies that <r(t) must remain finite. It follows 
from the above solutions that cr(t) can remain finite only if -y/3fi e > -jE^"> 

equivalently, 7 > ^7=5 - Under this condition, it can be checked easily that 

lim x ^ xe 3( - 1_7 ' )A = 0. Hence, f2 e = f2 as follows from ([191) . 

It can now be seen, for Q e = Q < |, that the quantities in (PD are all 
finite for tj < t < 00, implying that all the curvature invariants are finite. 
Thus, there is no singularity for tj < t < 00 when Q < §■ 

Consider the second case where K(t) in (|8]) may not remain < 1 for all 
t > tj. This is the case for the observed universe where 7 < and, hence, 
<r(t) grows as fast as or faster than t. Then there is a time, say t — t%, when 
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the value of K = 1. It follows that c < (cr(ti) + c) < and, hence, x(t) < 
for U < t < t\. Also, x(ti) is non vanishing since the case of x(£i) — > is 
same as the case described above. Therefore, we have that ^p(ii) = which 

implies, since 7^ 0> that ^(^i) = 0- This is a critical point of e A and is 
a maximum. 

For t > t\,K(t) becomes > 1 and, hence, A(t) < 0. Also, x < and 
Q(t) < 3. Then e = — 1 necessarily. Further evolution is analysed below. 

2 e. A(ti) < 0, < °> and < 3 

The above conditions imply that e = — 1, the constant c < in ([?]), and 
jfc{ti) > 0. Hence, K(ti) > 1. For i > ti, (cr(t) + c) which is negative at tj 
may or may not vanish for non zero X- 

Consider the case where (cr(t) + c) does not vanish and remains negative 
for non zero x- Therefore, x(£) < and, hence, x(t) decreases. For t > U, 
the scale factor e A decreases and, eventually, 7 can be taken to be > ~ when 
such "matter" is present. In fact, this is the case for the observed universe. 
Then, (1 — 37) < and a{t) decreases. Hence, (cr(t) + c) also decreases. 

It follows from (||) that K, initially > 1, is now decreasing since x an d 
e A are decreasing and (cr(t) + c) 2 is increasing. Its lowest value = ^ < |, 
achieveable at % = 0. Therefore, there exists a time, say t = tM, where 
K = 1 and, hence, ^ = 0. Clearly, x(^m) and x(£m) are non zero for the 

same reasons as in section 2c, implying that A(tjy) = 0. This is a critical 
point of e A and is a minimum. 

The existence of this critical point of e A can be seen in another way also. 
As (cr(t) + c) < grows in magnitude, it follows from equation (0) that 
fix 2 e~ 6A , whereas p is given by @. Note that x(t) < and, hence, x(t) 
decreases and eventually — > 0. Then, in the limit x ~^ 0> the term in 
equation @ dominates ^ term. Hence, in this limit where — > Q < | < 3, 

we have that A > 0. Since A < initially, this implies the existence of a zero 
of A. This is a critical point of e A and is a minimum. 

Note that all quantities remain finite for U < t < tj^- In particular, the 
quantities in ([LID are all finite, implying that all the curvature invariants are 
finite. Thus, there is no singularity for U < t < tj^- 

For t > t^j, we have A > 0, x < 0, < 3, and e = — 1. Further 
evolution then proceeds as described in section 2d. The evolution can thus 
become repetetive and oscillatory. 
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Consider now the case where (cr(t) + c) vanishes, say at t — t m /, for 
non zero X- Hence, x(t m ') vanishes and this is a minimum of X- F° r t > 
t m r, (cr(t) + c) and x are positive by continuity and, hence, x increases. Also, 
A < and the scale factor e A decreases. Eventually, 7 can be taken to be 
> I when such "matter" is present. In fact, this is the case for the observed 
universe. Then, (1 — 37) < and a(t) decreases. Hence, (cr(t) + c) also 
decreases. 

Now, as x increases and — > 1, (a(t) + c) may or may not vanish with 
X < 1- If (&(t) + c) vanishes with % < 1, then x vanishes. This is a critical 
point of x an d is a maximum, beyond which we have A < 0, x < 0, and 
e = — 1. Further evolution then proceeds as described in section 2c if Q > 3, 
and as in section 2e if f2 < 3. The evolution can thus become repetetive and 
oscillatory. 

If (<r(t)+c) does not vanish and remain positive for x < 1 then, eventually, 
X — > 1. It follows from equation (^) that fix 2 oc e _6A , whereas p is given by 
(|6p. Note that A < and e^ 4 is decreasing. Then, in this limit, the term 
in equation @ dominates or, if 7 = 1, is of the same order as the ^ term. 

Thus we have A < 0, x > 0, x — ^ 1> equivalently fi — > 00, and e = — 1. 
Also, -\ term in equation (f|) dominates or, if 7 = 1, is of the same order as 
the ^ term. But this is precisely the time reversed version of the evolution 
analysed in section 2b and in section 2a for the case which led to the initial 
conditions of section 2a, where the dynamics is clear in terms of and i/j{4>). 
Applying the results of sections 2a and 2b, it follows that <fi will cross the 
value after which and, hence, x begins to decrease. We then have A < 0, 
X < 0, and Q > 3. Also, e = — 1 necessarily. Further evolution then proceeds 
as described in section 2c. The evolution can thus become repetetive and 
oscillatory. 

3. Evolution of observed Universe 

We now use the results of the analysis in section 2 to describe the generic 
evolution of the observed universe. Note that the observed universe certainly 
contains dust (7 = 0) and radiation (7 = ~). As follows from (H), "matter" 
with larger 7 dominates the evolution for smaller value of e A and vice versa. 
Therefore, when e A is decreasing we take 7 > | eventually, and when e A 
is increasing we take 7 < eventually. We start with an initial time tj, 
corresponding to a temperature, say > 10 16 GeV, such that GUT symmetry 
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breaking, inflation, and other (matter) model dependent phenomena may 
occur for t > U only. Our observed universe is expanding at t{ and, hence, 
A(ti) > 0. For the sake of definiteness, we take x(U) > and Q(ti) > 3, 
equivalently ufa) > 0, as commonly assumed. 

We first describe the evolution for t > ti taking, as initial conditions 
A(tj) > 0, x(ti) > 0, and 0(t$) > 3. Then, e = 1 necessarily. To describe 
the evolution for t < ti, we reverse the direction of time and take, as initial 
conditions, A{ti) < 0, x(U) < 0, and VL(U) > 3. Then, e = — 1 necessarily. 
The required evolution is that for t > ti in terms of the reversed time variable, 
which is also denoted as t. The results of section 2 can then be applied 
directly. 



t > U 



Initially, we have A(U) > 0, x{U) > 0, and Q(ti) > 3. Then, e = 1 
necessarily. For t > ti, the evolution proceeds as described in section 2a. 
Both e A and x increase. Hence, fl increases. Eventually, as e A increases, we 
can take 7 < 0. Then, as t — > 00, e A and x evolve as given in (13). 

In particular, it can be seen that the quantities in (|TT| ) are all finite for 
ti < t < 00, implying that all the curvature invariants are finite. Hence, the 
evolution is singularity free. 

Thus, in the present day universe in our model, x(today) — > 1 and 
tt(today) -> 00. Also, as follows from (|), ±^(today) = -§|(1 - x) 4q_1 -> 
0. Therefore, our model satisfies the observational constraints imposed by so- 
lar system experiments, namely Q(today) > 2000 and -^^(today) < 0.0002. 

t < ti, equivalently (— t) > (—ti) 

To describe the evolution for t < t iy we reverse the direction of time. 
Then, in terms of the reversed time variable, also denoted as t, we have 
A{ti) < 0, x{ti) < O5 an d &{ti) > 3. Then, e = — 1 necessarily. For t > U, 
the evolution proceeds as described in section 2c. Both e A and x decrease. 
Hence, Q decreases. Then, as shown in section 2c, there exists a time, say 
t = t m > ti where K{t m ) = 1 in equation (||), e A ^ m) > 0, x(Pm) > 0, and 
X(t m ) < 0. Note that £l(t m ) < 3. 

As shown in section 2c, K{t m ) = 1 implies that A(t m ) = 0. Hence, the 
scale factor e A reaches a minimum. For t > t m , we then have A > 0, x < 0, 



12 



f2 < 3, and e = — 1. The evolution proceeds as described in section 2d. Now, 
however, the evolution for t > t m is complicated since the universe is known 
to contain "matter" with 7 = 0. Nevertheless, all qualitative features of its 
evolution can be obtained using the results of section 2. 

For t > t m , A > 0, x < 0; an d (c{t) + c) < 0. Hence, e A increases and 
X decreases and, eventually, "matter" with 7 < dominates the evolution. 
Then (1 — 37) > and a{t) begins to increase. Hence (cr(t) + c), which is 
negative, also begins to increase. For 7 < 0, a(t) grows atleast as fast as t, 
as follows from @. Then eventually, as described in section 2d, e A reaches 
a maximum eventually at time, say t = t\. Also, + c) < and, hence, 

x(ti) < 0. Note that this critical point of e A exists independent of the details 
of the model, as long as "matter" with 7 < is present. Further evolution 
then proceeds as described in section 2e. 

For t > ti, we then have A < 0, x < 0, and Vt < 3. Also, (a(t) + c) < 
but a(t) remains increasing. Hence, depending now on the details of the 
model, (er(t) + c) may remain negative for all t > t\, or it may reach a zero 
and become positive. 

In the first case, (cr(t) + c) remains negative for all t > t\. Hence, x < 
and x continues to decrease. The scale factor e A , as shown in section 2e, 
reaches a minimum at time, say t = tj\t > t\, with x{^m) > non vanishing. 
For t > Imi e A increases, and we have A > 0, x < 0, and < 3. Further 
evolution then proceeds as described in section 2d. 

Note that this means that the scale factor increases and reaches a maxi- 
mum, then decreases and reaches a minimum, then increases and so on. The 
existence of maxima of e A is model independent as long as "matter" with 
7 < is present, which is the case for the observed universe. The minima of 
e A are all non zero for the reasons described in section 2c. Their existence de- 
pends on whether (a(t) + c) remains negative or not, but is otherwise model 
independent as long as "matter" with 7 > | is present, which is the case 
for the observed universe. The evolution can thus become repetetive and 
oscillatory. 

In the second case, (<r(t) + c) reaches a zero at time, say t = t m > > t\ and 
becomes positive. It then follows that x(t m ') — and that this is a minimum 
of x- For t > t m /, we have A < 0, x > 0, and f2 < 3. Hence, e A decreases 
and x increases. Eventually, we can take 7 > |. 

Now, (1 — 37) < and a(t) begins to decrease or remains constant. Hence 
(cr(t) + c), which is positive, also begins to decrease or remains constant. 
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Thus, depending on the details of po an d 7, {&{t) + c) may or may not vanish 
with x < 1- If i a {t) + c ) vanishes with x < 1, then % vanishes. This critical 
point is a maximum of x, beyond which we have A < 0, x < 0, and e = — 1. 
Further evolution then proceeds as described in section 2c if Q > 3 and as in 
section 2eifQ < 3. The evolution can thus become repetetive and oscillatory. 

If (a(t)+c) does not vanish and remain positive for x < 1 then, eventually, 
X — > 1. It follows from equation (^) that f2x 2 oc e~ 6A , whereas p is given by 
@. Note that A < and e A is decreasing. Then, in this limit, the term 
in equation (Q) dominates or, if 7 = 1, is of the order of the ^ term. 

Thus we have A < 0, x > 0) X 1> equivalently — > 00, and e = — 1. 
Also, ^p- term in equation (|4]) dominates or is of the order of the ^ term. 
But this is precisely the time reversed version of the evolution analysed in 
section 2b, where the dynamics is clear in terms of and ip((f>). Applying the 
results of section 2b, it then follows that (f> will cross the value after which 

and, hence, x w iU begin to decrease. We then have A < 0, x < 0, f2 > 3. 
The evolution then proceeds as described in section 2c. The evolution can 
thus become repetetive and oscillatory. 

Thus it is clear that depending on the details of the model, the universe 
undergoes oscillations, perhaps infinitely many. As follows from the above 
description, the oscillations can stop, if at all, only in the limit x — > 0. The 
solutions then are given by (|2l|) - (p3|) . In particular, however, the quantities 
in ( |TT| ) all remain finite during the oscillations, implying that all the curvature 



invariants remain finite. They also remain finite in the solutions (^l|) - ( p3|) 
ifO o <|in(0). 

Thus, it follows that if < | then the quantities in (p7|) are all finite 
for ti < t < 00, implying that all the curvature invariants are finite. Hence, 
the evolution is singularity free. 

In summary, the evolution of a realistic universe, such as our observed 
one, with the initial conditions A(tj) > 0, > 0, and f2(tj) > 3 proceeds 
in the present model as follows. For t > ti, the scale factor e A increases 
continuously to 00. The field x increases continuously to 1. Correspondingly, 
Q increases continuously to 00. 

For t < ti, the scale factor e A decreases and reaches a non zero mini- 
mum. It then increases and reaches a maximum, then decreases and reaches 
a minimum, then increases and so on, perhaps ad infintum. The oscillations 
can stop, if at all, only in the limit x ~~ * 0. The solutions then are given by 
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( |2TD - (p3|). The field x ma Y; depending on the details of model, continuously 
decrease to 0, or undergo oscillations, its maxima always being < 1. 

Also, the curvature invariants all remain finite for — oo < t < oo. Hence, 
the evolution, although complicated and model dependent in details, is com- 
pletely singularity free if fi(x) satisfies @. Thus we have that a homogeneous 
isotropic universe, such as our observed one, evolves with no big bang or any 
other singularity in a class of models where fi(x) satisfies @). The time 
continues indefinitely into the past and the future, without encountering any 
singularity. 

When the initial conditions are different, the evolution can again be anal- 
ysed along similar lines. However, the main result that the evolution is 
singularity free remains unchanged. 
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